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摘要 :利用 黄 淮海 平原 45 个 气象 站 点 1961 一 2014 年 月 值 气象 数据 ,基于 Penman-Monteith 蒸 散 模 型 计算 了 标准 化 降水 蒸 散 指 
数 (SPEI) ,对 黄 淮海 平原 近 34 年 干旱 变化 趋势 发生 频 率 和 持续 性 特征 进行 了 分 析 ,并 探讨 了 SPEI 指数 与 河南 .河北 和 山东 省 
农业 干旱 面积 的 关系 ,结果 表明 :(1) 改 用 Penman-Monteith 芋 散 公式 后 ,SPEI 干旱 指数 在 黄 淮 海平 原 呈 整体 上 升 趋势 , 即 趋 于 
湿润 ;(2) 近 54 年 干旱 演变 具有 明显 的 年 代 际 差异 ,20 世纪 60 年 代 干 旱 频率 最 高 ,而 21 世纪 初 (2000 一 2014) 干旱 频率 整体 候 
低 ;(3) 黄 淮海 平原 干旱 发 生 具 有 持续 性 的 特点 ,20 世纪 60 年 代 遭 受 的 持续 性 干旱 最 为 严重 ,平均 干旱 持续 时 长 约 2.6 个 月 ,21 
世纪 初 下 降 到 1.5 个 月 ;(4) 河 南 .河北 和 山东 省 的 农业 干旱 面积 年 际 变化 表明 ,干旱 面积 呈 减 少 趋势 ,2000 年 以 后 年 均 受 灾 面 
积 .成 灾 面 积 和 绝收 面积 比 2000 年 之 前 分 别 下 降 了 58.0% .44.4% 和 49.1%;(5) 农 业 干 旱 面积 与 SPEI 具有 中 等 以 上 的 相关 强 
度 , 其 中 对 山东 省 受灾 成 灾 和 绝收 面积 相关 系数 + 达到 -0.7 以 上 ,表明 基于 Penman-Monteith 蒸 散 模 型 的 SPEI 指数 在 黄 淮 海 
平原 具有 良好 的 适用 性 。 
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Abstract: The quantitative analysis of drought is crucial essential for drought risk assessment. The standardized 
precipitation evapotranspiration index (SPEI) has been widely used as an effective approach to quantitatively analyze the 
trends, duration, frequency, and severity of drought. However, how to calculate the evapotranspiration is big challenge for 
the reliable and accuracy of SPEI results. The most previous studies calculated SPEI based on an empirical 
evapotranspiration equation, such as the Thornthwaite method, rather than the physical Penman-Monteith equation. 
Moreover, majority of the studies analyzed the spatio-temporal pattern of the index without linkage of the relationship 
between the climactic drought index and actual cropland drought areas. This study utilized SPEI based on the Penman- 
Monteith equation to explore the variations of drought frequency and duration during 1963—2014 in 3H Plain, and 


established the relationship between SPEI and the actual field drought areas by Pearson correlation in Henan, Hebei, and 
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Shandong Province separately, which aims to understand the trends of drought during the past 50 years and to explore the 
relationship between the climatic drought index and actual field drought areas. The results showed that the SPEI had upward 
trend for 1, 3, 6, and 12-month scales in most areas of the 3H Plain while as the 1 -month and 3-month scales, the wetting 
trend was significant ( P<0.1) in the northern regions of 3H Plain. Drought frequency in the 1960s was the highest during 
the past 54 years, but lowest in 2000 一 2014. The longest duration was occurred in the 1960s by comparing the decadal 
spatial distribution of drought duration. Moreover, the mean drought duration declined from 2.6 months in the 1960s to 1.5 
months during 2000 一 2014. The changes of annual drought areas decreased in Henan, Hebei, and Shandong Province, 
which could be partly attributed to the lower frequency and shorter duration of drought detected by SPEI. Correlation 
analyses indicated that the index series of Dec-SPEI- 12 had medium to high correlation with observed cropland drought 
areas. For example, Pearson’ s r between Dec-SPEI- 12 and the drought area of Shandong Province were -0.7，-0.7，and 
-0.8 for affected, disaster and no harvest drought areas, respectively. The results of the study suggest a wetting trend in the 
3H Plain during the past 54 years, and the chosen of Penman-Monteith evapotranspiration model with decreasing 
evapotranspiration over the past years in the 3H Plain attributes to this wetting trend. The high correlation between the 
climatic drought index and actual drought area indicated that SPEI could be used as a reference or trigger for establishing a 


drought warning system in the 3H Plain. 
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干旱 灾害 影响 范围 广 持 续 时 间 长 ,一旦 发 生 所 造成 的 影响 大 ,因而 备 受 关注 ,是 当前 气候 变化 研究 的 热 
点 领域 '"。 我 国 地 处 东亚 季风 区 ,地 理 条 件 复杂 ,气象 灾害 频 发 '"| ,而 旱灾 是 我 国 当 前 最 主要 的 气象 灾害 之 
一 ,对 农业 生产 .水 资源 .生态 环境 和 人 类 活动 产生 重要 影响 。1961 一 2010 年 ,我 国 干旱 年 均 受 灾 面 积 约 
23230 千 hm? , 约 占 农 作物 播种 总 面积 的 15.6%; 其 中 成 灾 面 积 约 10953 千 hm? , 约 占 播种 总 面积 的 7.3%D1 。 
IPCC 报告 以 及 我 国 气候 变化 预 估 分 析 均 表明 ,气候 变化 将 会 造成 极端 气候 事件 及 气候 灾害 的 频率 和 强度 明 
显 增强 ,气候 灾害 影响 的 区 域 增多 ,进而 对 社会 经 济 可 持续 发 展 造成 不 利 影响 。 因 此 定量 化 研究 干旱 演变 特 
征 对 制定 适应 气候 变化 策略 .干旱 风险 管理 具有 科学 指导 意义 。 

20 世纪 以 来 ,针对 干旱 影响 的 部 门 和 地 区 ,特殊 的 应 用 和 对 其 理解 的 程度 不 同 ,国际 上 提出 了 许多 干旱 
定量 化 研究 方法 。 其 中 ,基于 复杂 水 分 平衡 模式 提出 的 Palmer 干旱 指数 (PDSI) 是 干旱 研究 发 展 史 上 的 里 程 
碑 , 有 效 地 解决 了 区 域 干旱 监测 、 预 报 难题 ,但 是 缺乏 时 空 比较 的 有 效 性 “* 引 。McKee 等 人 !9 提出 标准 化 降水 

旨 数 (SPI) 来 蔡 代 Palmer 指数 ,SPI 具有 时 间 尺 度 灵 活 的 特点 ,计算 简便 ,但 仅 考虑 降水 因素 ,忽略 了 水 分 平 
衡 的 影响 ,对 研究 气候 变 暧 背景 下 干旱 发 生 趋 势 效 果 不 佳 '" 。 基 于 此 ,2010 年 Vicente-Serrano 等 人 [1 提出 了 
标准 化 降水 莹 散 指 数 (SPEI) ,该 指数 基于 降水 和 莹 散 的 差额 , 既 保 留 了 PDSI 对 温度 敏感 的 特点 ,又 具备 SPI 
适合 多 尺度 、. 多 空间 比较 的 优点 , 适 于 全 球 变 暖 背景 下 干旱 特征 的 分 析 '5。2014 年 , Beguerfa S 与 Vicente- 
Serrano 等 人 对 SPEI 指数 进行 了 回顾 '” ,指出 由 于 不 同 落 散 公式 物理 意义 的 不 同 、 考 虑 的 气象 要 素 不 同 ,计算 
的 SPEI 指数 会 产生 差异 。Thornthwaite 蒸 散 模型 所 需 气象 要 素 仅 有 温度 ,计算 简便 '” ,基于 该 模型 的 SPEI 得 
到 广泛 应 用 "2 。 目 前 采用 基于 Penman-Monteith 模型 的 SPEI(SPEI-PM) 对 干旱 的 研究 较 少 , Ming Bo 等 
人 ' 引 分 析 了 SPEI-PM 指数 与 华北 平原 夏 玉米 产量 的 关系 ;Kai Xu 等 人 "利用 SPEI-PM 指数 对 我 国 1961 一 
2012 年 的 干旱 特征 进行 分 析 。 因 此 , 改 用 应 用 范围 更 广 物 理 意 义 更 明确 的 Penman-Monteith 模型 后 干旱 特 
征 值 得 探讨 。 

黄 淮海 平原 是 我 国 气候 变化 敏感 区 与 脆弱 区 ,气候 变 暧 导致 降水 变异 程度 加 大 ,加 之 地 表 水 时 空 分 布 不 
均 .地 下 水 超 采 严重 55) ,干旱 风险 趋向 不 稳定 。 对 该 区 干旱 的 研究 成 果 丰 富 ,很 好 地 揭示 了 干旱 的 南北 差异 
和 季节 性 特点 ,但 许多 研究 仅 对 干旱 指标 进行 时 空 分 析 , 而 气象 干旱 指标 与 实际 农业 旱情 的 联系 不 清 。 因 此 ， 
本 文 改 用 基于 Penman-Monteith 模型 的 SPEI 指数 对 黄 淮 海平 原 1961 一 2014 年 的 干旱 发 生 特 征 进行 分 析 , 从 
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干旱 指数 时 空 特征 和 演变 趋势 出 发 ,寻求 干旱 灾害 变化 规律 ;从 历史 灾情 数据 与 干旱 指数 的 关系 出 发 ,验证 气 
象 干旱 指数 与 农业 实际 干旱 的 关系 ,为 该 地 区 灾害 风险 特征 和 灾害 风险 评估 提供 依据 。 


1 材料 与 方法 


1.1 区 域 概况 

黄 淮 海平 原 位 于 燕山 以 南 .淮河 以 北 (112°33' 一 120°17 祷 ,31°14' 一 50°25'N) ,由 黄河 ,淮河 和 海河 冲积 平 
原 及 部 分 丘陵 山区 组 成 , 属 半 湿润 地 区 ,年 降水 量 500 一 800 mm, 呈 南 多 北 少 的 分 布 格局 ,年 潜在 蒸 散 量 为 
1000 mm 左右 ,大 部 分 区 域 降水 处 于 亏 缺 状态 ,是 我 国 的 干旱 重 灾区 "之 一 。 同 时 , 黄 淮 海平 原 是 我 国 的 重 
要 粮食 生产 基地 ,主要 种 植 方 式 是 冬小麦 - 夏 玉米 轮作 ,小 麦 和 玉米 种 植 面积 分 别 占 全 国 的 61% 和 31%。 本 文 
根据 中 国 农 作 制 ' 趾 黄 淮 海平 原 半 湿 润 暧 温 灌 溉 集约 农 作 区 划 , 将 黄 淮 海平 原 分 为 以 下 6 个 农业 亚 区 :工区 燕 
山 太行 山 山 前 平原 水 浇 地 二 熟 区 ; [区 环 渤 海滨 海外 向 型 二 熟 农 汐 区 ; 耳 区 海河 低 平原 缺 水 水 浇 地 二 熟 兼 旱 
地 一 熟 区 ;区 鲁 西 平原 水 浇 地 二 熟 兼 一 熟 区 ; V 区 黄 淮 平原 南阳 伟 地 水 浇 地 旱地 二 熟 区 ; VI 区 江淮 平原 麦 稳 
二 熟 区 ,如 图 1。 
1.2 ”资料 来 源 

研究 区 内 45 个 站 点 的 1961 一 2013 年 月 值 气象 数据 来 自 国 家 气象 局 ,包括 降水 (mm) .最 高 气温 (% ) .最 
低 气 温 (% ) .风速 (m/s) 日照 时 数 (h) 、 相 对 湿度 (%) 经 纬度 和 海拔 高 度 。 研 究 区 及 所 选 气象 站 点 分 布 见 
图 1。 历 史 干 旱灾 情 数据 来 源 于 农业 部 种 植 业 管理 司 灾情 数据 库 。 本 文 收 集 了 时 间 序 列 较 长 的 河北 河南 以 
及 山东 省 历年 干旱 受灾 成 灾 和 绝收 面积 数据 ,其 中 受灾 面积 序列 为 1971 一 2013 年 ,成 灾 面 积 序列 为 1971 一 
2012 年 ,绝收 面积 为 1982 一 2012 年 。 


114° 116° 118° 120° 122° E 


图 1 黄 淮海 平原 气象 站 点 分 布 图 


Fig.l The location of meteorology stations in Huang-Huai-Hai Plain 
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1.3 SPFI 指数 
SPEI 指数 通过 计算 月 降水 与 潜在 蒸 散 量 的 差 值 并 进行 正 态 标准 化 处 理 得 到 。 
第 一 步 ,计算 逐 月 潜在 蒸 散 量 BTu。 本 文采 用 FAO-56(1998 ) 推荐 的 Penman-Monteith 模型 计算 : 


900 
0.408A( 民 ， C) 十 Yy pr ey) 


A+y(l1 +0.34U,) Wy 
式 中 ,£7, 为 潜在 燕 散 量 (mm) ; 4 为 温度 随 饱 和 水 汽 压 变 化 的 斜率 ( kPa/%C ) ; U, 为 离 地 2 m 高 处 风速 (my 
s) ; 6, 为 空气 饱和 水 汽 压 (kPa) ; e, 为 空气 实际 水 汽 压 (kPa) ; 7 为 平均 气温 (% ) ; y 为 湿度 表 常 数 ( kPa/ ) ; 
RR, 为 到 达 地 面 的 净 辐 射 (MJ m”d”) ; G 为 土壤 热 通 量 密度 (MJ md ) 。 
第 二 步 , 计 算 逐 月 降水 量 与 蒸 散 量 的 差额 : 

D,=P,- (ET,), (2) 
式 中 , D; 是 计算 的 时 间 尺 度 内 降水 与 蒸 散 差额 的 累计 值 , Di 为 第 i 年 第 j 个 月 开始 ,下 个 月 内 的 累积 降水 蒸 
散 差额 。 


ET, = 


Di = Ee J < 天 
p= FD, 了 三 大 
第 三 步 ,对 D; 数据 序列 进行 拟 合 。 研 究 发 现 '" ,采用 三 参数 的 Log-Logistic 概率 分 布 函数 的 拟 合 效果 
最 好 。 
4 一 P 1 % 一 pl 
N= L(t) Ex ”| (3) 
[0.4 CQ (64 
式 中 ,参数 a ,8 , y 可 以 采用 线性 和 矩 (L-moment ) 方 法 拟 合 获得 : 
2w1 — wo (oo - 2w;)B _ rrl+1 rr1_1 
CE 
T(8) 是 关于 B 的 Gamma 函数 。 由 此 可 以 得 到 D, 的 概率 密度 的 累计 概率 密度 函数 : 
B1] -1 
Feo)= 1+ ~ ) | (4) 
x%—y 
第 四 步 , 对 累计 概率 密度 进行 正 态 标准 化 。 超 过 某 个 D, 值 的 概率 为 P = 1 - F(x) ,概率 加 权 和 矩 w = 
V -2In(P) 。 


当 累 计 概 率 P <0.5 时 ,采用 (5) 式 计算 SPEI。 
Co 十 CIO 十 CO 
SPEI = w 一 5 (5) 
l+diw + dw + dw 
式 中 , d=1.432788; d, =0.189269; d; =0.001308; co=2.515517; c, =0.802853; c, =0.013028。 


当 P >0.5 时 ,采用 (6) 式 计算 SPEI。 


C0 tc + co 
SPET =- (ow- a ;) (6) 
l+dw + dw + dw 


SPEI 指数 是 具有 多 时 间 尺 度 的 标准 化 干旱 指标 ,等 级 划分 及 对 应 的 发 生 概率 '" 见 表 1。 不 同时 间 尺 度 
的 SPEI 反映 不 同 水 源 的 可 利用 性 ,土壤 湿度 对 短 时 间 尺 度 SPEI 较为 敏感 ,地 下 水 河川 径流 和 水 库 划 水 量 对 
长 时 间 尺 度 SPEI 较为 敏感 (中 。 本 文 利用 SPEI-1、SPEI-3、SPEI-6 和 SPEI-12 分 别 表示 较 短 时 间 (1 个 月 .3 个 
月 .6 个 月 ) 和 较 长 时 间 (12 个 月 ) 的 时 间 尺 度 SPEI 研究 干旱 趋势 随时 间 的 变化 规律 。 
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1.4 数据 处 理 

本 文采 用 Beguerfa S 等 开发 的 R 语言 扩展 包 计 算 SPEI ( https://cran. r-project. org/ web/packages/ 
SPEL)'” ,并 选取 PM 模型 为 潜在 蒸 散 模型 ,计算 了 研究 区 45 个 气象 站 点 1963 一 2014 年 逐 月 多 尺度 SPEI 
值 , 对 黄 淮海 平原 的 干旱 趋势 随时 间 变 化 进行 分 析 : 


表 1 SPEI 干旱 等 级 划分 
Table 1 Drought classifications based on SPEI 


干旱 等 级 概率 /% SPEI 值 干旱 等 级 概率 /9% SPEI 值 
Drought class Probability SPEI value Drought class Probability SPEI value 
极端 湿润 Extremely wet 2.3 =2.0 中 度 干 旱 Moderately dry 9.2 -1.5 一 -1.0 
严重 湿润 Severely wet 4.4 1.5 一 2.0 严重 干旱 Severely dry 4.4 -2.0 一 -1.5 
中 度 湿润 Moderately wet 9.2 1.0 一 1.5 极端 干旱 Extremely dry 2.3 二 -2.0 
正常 年 份 Near normal 68.2 -1.0 一 1.0 


二 干旱 指数 的 趋势 分 析 : 采 用 Mann-Kendall( MK ) 趋势 检验 方法 对 黄 淮海 平原 45 个 气象 站 点 年 均 SPEI-1、 
> SPEI-3、SPEI-6、SPEI-12 的 变化 趋势 显著 性 进行 检验 。MK 趋势 检验 是 非 参 数 统计 检验 方法 ,计算 简单 ,是 目 
前 比较 常用 的 趋势 诊断 方法 ,计算 方法 参考 文献 .| 。 
干旱 频率 :按照 表 1 的 干旱 等 级 划分 方法 对 黄 淮海 平原 6 个 农业 亚 区 1963 年 1 月 一 2014 年 12 月 的 
SPEI-1 SPEI-3 SPEI-6、 SPEI-12 的 干旱 频率 p 进行 计算 ,方法 如 下 : 
p=n/N x 100% (7) 
污 式 中 ,p 表示 干旱 发 生 频 率 ; n 为 数据 序列 中 干旱 发 生 的 次 数 ; N 为 数据 序列 数 。 
a 干旱 的 持续 性 分 析 :根据 表 1 对 干旱 等 级 的 划分 ,本 文 定义 SPEI-3 大 -1.0 的 连续 月 份 为 持续 时 长 。 
OO SPEI 与 干旱 面积 的 分 析 :为 探讨 气象 干旱 指标 与 实际 农业 干旱 的 关系 ,采用 Pearson 相关 系数 法 分 析 灾 
“” 情 面积 与 各 时 间 尺 度 SPEI 指数 的 关系 。Potopova V 等 人 使 用 该 方法 探讨 了 捷克 共和 国 不 同 农 作物 气候 产量 
与 SPEI 的 相关 关系 。 


2 结果 与 分 析 


志 2.1 干旱 指数 的 趋势 变化 
© 图 2 描述 了 黄 淮海 平原 干旱 指数 的 变化 趋势 ,其 中 MK 趋势 系数 小 于 0 表明 指数 呈 下 降 趋势 , 即 干旱 化 ， 
反之 , 则 趋 于 湿润 ;并 且 趋 势 系数 的 绝对 值 大 于 1.64(1.96) 时 ,表明 变化 趋势 的 通过 90% (95%) 置信 水 平 。 
从 趋势 系数 的 分 布 可 以 看 出 ,基于 PM 蒸 散 模型 的 SPEI 在 各 时 间 尺 度 上 表现 为 上 升 趋势 站 点 数 大 于 下 降 趋 
势 的 站 点 数 , 即 黄 淮海 平原 整体 呈 湿 润 趋势 。 年 均 SPEI-1 上 升 趋势 最 为 明显 (图 2) ,有 13 个 站 点 达到 了 
90% 以 上 显著 性 水 平 ,其 中 11 个 站 点 达到 了 95% 置 信 水 平 , 仅 有 3 个 站 点 呈 微 弱 的 下 降 趋势 ;年 均 SPEI-3 有 
9 个 站 点 达到 了 90% 以 上 显著 性 水 平 (图 2) ,其 中 7 个 站 点 达到 了 95% 和 置信 水 平 , 仅 有 7 个 站 点 呈 微 弱 下 降 
趋势 ;从 SPEI-6 和 SPEI- 12 来 看 (图 2) , 黄 淮海 地 区 干旱 指数 呈 上 升 趋势 的 站 点 分 别 占 75.6% 和 73.3% ,其 中 
各 有 6 和 7 个 站 点 达到 90% 以 上 显著 性 水 平 , 分 别 有 11 和 12 个 站 点 呈 微 弱 下 降 趋 势 。 从 变化 趋势 的 空间 差 
异 来 看 ,湿润 化 站 点 主要 集中 在 黄 淮海 平原 中 北部 ,例如 河北 省 廊坊 沧州 ,秦皇岛 、 乐 亭 和 南宫 站 的 SPEI-1 
和 SPEI-3 指数 上 升 趋势 均 达 到 95 多 的 显著 性 ,山东 与 河南 交界 的 朝阳 、 充 州 . 尚 丘 和 西 华 站 的 SPEI- 1 和 
SPEI-3 趋势 同样 都 达到 了 显著 性 ;干旱 化 站 点 分 布 较为 零散 ,主要 集中 在 京 津 和 崇 北 地 区 ,但 并 未 达到 显 
著 性 。 
2.2 干旱 演变 特征 

为 了 分 析 黄 淮海 平原 1963 一 2014 年 干 湿 演变 特征 , 本文 利用 气象 中 常用 的 哈 默 图 方法 ( Hovmoller-type 
diagrams) 对 黄 淮海 平原 6 个 农业 亚 区 多 时 间 尺 度 的 干旱 演变 进行 分 析 ( 图 3) ,能 够 直观 地 显示 干旱 发 生 的 起 
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图 2 黄 淮海 平原 45 个 气象 站 点 年 均 SPEI-1、SPEI-3、SPEI-6 和 SPEI-12 变化 趋势 


Fig.2 Trend variations of annual SPEI with time scales of one-month three-month, six-month and twelve-month in 3H Plain 


止 时 间 ,严重 程度 和 持续 时 长 。 

从 图 3 可 以 看 出 , 短 时 间 尺 度 ( 如 SPEI-3) 的 干 湿 指数 波动 频繁 ,反映 了 短期 降水 对 干旱 程度 的 影响 ; 随 
着 时 间 尺 度 加 大 ,波动 周期 相对 较 长 ,体现 了 干 湿 的 季节 性 变化 规律 ;长 时 间 尺 度 ( 如 SPEI- 12) 干 湿 变 化 更 为 
稳定 ,波动 周期 长 ,更 能 反映 干旱 的 年 际 变化 特征 。 从 干旱 指数 的 演变 特征 来 看 , 呐 淮 海平 原 6 个 农业 亚 区 干 
旱 发 生 具 有 明显 的 年 代 际 特征 。 对 于 工 一 焉 区 (图 3) ,干旱 主要 发 生 在 20 世纪 60 年 代 中 后 期 (1965 一 


1968) .20 世纪 


80 年 代 初 期 (1980 一 1984) .20 世纪 90 年 代 末 (1998 一 1999) 和 21 世纪 2002 年 左右 ;IV 区 (图 


3) 有 5 条 明显 的 干旱 带 ,分 别 是 1966 一 1969 .1981 一 1982 .1989 一 1991 .1996 一 1999 和 2002 一 2003;V 和 WI 区 
(图 3) 具 有 相似 的 干旱 分 布 ,为 1966 一 1968 .1978 一 1980 .1995 一 1996 .2001 一 2002 .2010 一 2011 和 2014 年 左 
右 。 这 与 我 国 历史 干旱 资料 记录 较为 一 致 , 荣 艳 淑 等 ”的 研究 结果 指出 的 典型 干旱 年 份 如 1965 年 .1972 年 、 
1986 年 .1997 年 及 2001 年 ,以 及 几 个 典型 干旱 时 期 如 1965 一 1967 年 .1980 一 1981 年 .1991 一 1992 年 .1999 一 
2002 年 以 及 2006 一 2007 年 在 SPEI 指数 时 间 序 列 中 均 得 到 较 好 的 体现 ,反映 了 SPEI 指数 在 华北 地 区 旱 洲 趋 
势 分 析 中 具有 较 好 的 适用 性 。 
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图 3 黄 淮 海平 原 各 农业 亚 区 多 时 间 尺 度 干 旱 历 史 演 变 
Fig.3 Spatiotemporal evolution of the SPEI series from 1 to 24 month scales at different sub-region 


为 定量 揭示 黄 淮 海平 原 各 农业 亚 区 干旱 发 生 的 年 代 差 异 , 本 文 计算 了 各 亚 区 不 同年 代 干 旱 (SPEI<-1) 发 
生 的 频率 , 见 表 2。 总 体 上 看 , 黄 淮 海平 原 干旱 频率 最 高 的 年 代 是 1960s,21 地 纪 以 后 干旱 频率 下 降 明 显 。 
区 干旱 频率 最 高 的 年 代为 1960s 和 1970s,4 个 时 间 尺 度 下 1960s 干旱 频率 分 别 为 23.00% 、21.88%、16.67% 和 
21.88% ,而 2000 年 以 后 的 干旱 频率 最 低 , 分 别 为 13.33% 、10.00% .8.89% 和 7.22%; 工区 干旱 频率 最 大 的 年 代 
仍旧 是 1960s ,各 时 间 尺 度 干 旱 频率 分 别 为 18.75% 、23.96%、16.67% 和 14.58%,2000 年 以 后 最 低 , 分别 为 
8.89% .9.44% .9.44% 和 7.22%; 焉 区 1960s 各 时 间 尺 度 的 干旱 频率 明显 高 于 其 他 年 代 ,依次 为 27.08% 、 
26.04% 28.13% 和 29.17% ,2000 年 以 后 干旱 频率 仍然 最 低 ; 区区 1960s 和 1980s 干旱 频率 差异 不 大 ,1970s 和 
2000 年 以 后 干旱 频率 最 低 ; V 区 1960s 干旱 频率 依然 明显 高 于 其 他 年 代 , 不 同时 间 尺 度 干旱 频率 分 别 为 
26.04% .26.04% .20.83% 和 26.04% ,2000 年 以 后 干旱 频率 依旧 处 于 最 低 ;对 于 VI 区 ,1960s 和 1970s 干旱 频率 
差异 不 大 ,但 是 1980s 干旱 频率 明显 低 于 其 他 年 代 ,分 别 只 有 11.67% .12.50% .9.17% 和 3.33% 。 
2.3 干旱 持续 性 特征 

对 黄 淮 海平 原 45 个 气象 站 点 干旱 最 长 持续 时 间 的 年 代 际 差异 进行 分 析 ,结果 表明 , 黄 淮海 平原 干旱 持续 
性 具有 明显 的 年 代 际 特征 (图 4) 。1960s 年 代 是 遭受 持续 性 干旱 影响 最 为 严重 的 10 年 ,最 长 干旱 持续 期 达到 
6 个 月 以 上 的 站 点 占 57.8% ,其 中 有 4 个 站 点 最 长 持续 时 间 达 到 9 个 月 ,分 别 是 天 津 、 沧 州 . 南 宫 和 黄骅 ,持续 
的 时 段 缘 为 1968 年 2 月 一 1968 年 10 月 (图 4) ;而 1970s 和 1980s 持续 时 间 达 到 6 个 月 以 上 的 站 点 数 分 别 为 
11 和 8 个 , 且 各 有 1 个 站 点 持续 时 长 达到 11 个 月 ,分 别 为 廊坊 (1975 年 3 月 一 1976 年 1 月 ) 和 塘 洁 (1988 年 
12 月 一 1989 年 10 月 ) (图 4) ;90 年 代 仅 有 1 个 站 点 达到 6 个 月 持续 长 度 (图 4) ,2000 一 2014 年 这 15 a 间 , 也 
仅 有 6 站 点 达到 了 6 个 月 持续 时 长 (图 4) 。 黄 淮海 平原 45 个 气象 站 点 干旱 平均 持续 月 数 的 年 际 变化 呈 波 动 
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下 降 趋势 (图 4) ,1960s 的 平均 持续 时 长 最 长 ,达到 2.6 个 月 每 次 ;1970s 和 1980s 逐步 下 降 ,分 别 为 1.9 和 1.7 
个 月 每 次 ;1990s 有 所 回升 ,为 2.0 个 月 每 次 ;21 世纪 初 最 低 (2000 一 2014) ,为 1.5 个 月 每 次 。 因 此 ,从 干旱 最 
长 持续 时 长 和 平均 持续 时 长 的 年 代 际 变化 来 看 ,1960s 是 遭受 持续 性 干旱 最 为 严重 的 10 年 ,而 本 世纪 初 干旱 
持续 性 有 所 减弱 。 


表 2， 黄 淮海 平原 各 亚 区 干旱 频率 的 年 代 际 变化 特征 /% 


Table 2 Decadal Variation of drought frequency at different scales for each region in 3H Plain 


亚 区 /Sub-regions 1960s 1970s 1980s 1990s 2000 一 2014 
工区 /Sub-region I 
SPEI-1 25.00 24.17 15.00 18.33 13.33 
SPEI-3 21.88 18.33 8.33 11.67 10.00 
SPEI-6 16.67 16.67 13.33 11.67 8.89 
SPEI- 12 21.88 21.67 16.67 16.67 7.22 
工区 /Sub-region I[ 
= SPEI- 1 18.75 9.17 10.83 14.17 8.89 
pe SPEI-3 23.96 11.67 11.67 14.17 9.44 
SPEI-6 16.67 5.83 10.83 14.17 9.44 
SPEI- 12 14.58 2.50 10.83 17.50 7.22 
亚 区 /Sub-region 亚 
SPEI-1 27.08 22.50 13.33 13.33 15.56 
SPEI-3 26.04 20.00 16.67 14.17 10.56 
SPEI-6 28.13 19.17 15.83 15.00 7.22 
SPEI- 12 29.17 10.00 8.33 19.17 8.33 
IV 区 /Sub-regionJV 
SPEI-1 20.83 19.17 17.50 20.83 15.00 
SPEI-3 19.79 13.33 20.00 17.50 12.22 
SPEI-6 17.71 11.67 27.50 19.17 13.89 
SPEI- 12 31.25 4.17 35.83 20.83 12.78 
之 V 区 /Sub-region V 
一 SPEI- 1 26.04 15.00 12.50 16.67 14.44 
SPEI-3 26.04 11.67 13.33 14.17 11.67 
SPEI-6 20.83 11.67 14.17 12.50 10.56 
SPEI- 12 26.04 11.67 18.33 11.67 9.44 
VW 区 /Sub-region VI 
SPEI- 1 21.88 20.00 11.67 18.33 13.89 
SPEI-3 27.08 21.67 12.50 18.33 15.00 
SPEI-6 18.75 21.67 9.17 19.17 12.22 
SPEI- 12 26.04 19.17 3.33 22.50 14.44 


对 黄 淮 海平 原 所 有 干旱 事件 进行 排序 ,分 析 持 续 时 间 最 长 的 10 次 干旱 事件 , 见 表 3。 廊 坊 和 塘 洛 分 别 在 
1975 年 3 月 一 1976 年 1 月 和 1988 年 12 月 一 1989 年 10 月 发 生 了 持续 11 个 月 的 干旱 事件 , 且 都 达到 了 中 旱 
以 上 水 平 (SPEI<-1.5) ;另外 ,从 1968 年 2 月 至 1968 年 10 月 ,河北 中 北部 的 石家庄 、 黄 骅 、 南 宫 和 沧州 以 及 天 
津 同 时 发 生 了 持续 9 个 月 左右 的 中 到 大 旱 ,表明 1968 年 华北 地 区 特别 是 京 津 惠 地 区 发 生 了 一 次 严重 的 干旱 
事件 。 

2.4 SPEI 指数 与 干旱 受灾 成 灾 和 绝收 面积 的 关系 

黄 淮海 平原 是 我 国 粮食 的 主 产 区 ,干旱 的 发 生 直接 对 该 地 区 的 农业 生产 造成 影响 。 为 了 探讨 SPEI 指数 

与 农业 干旱 的 关系 ,本文 收集 了 资料 连续 性 较 好 .序列 较 长 的 河北 .河南 以 及 山东 3 省 历年 干旱 受灾 .成 灾 和 
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图 4 黄 淮海 平原 各 年 代 干旱 最 长 持续 月 数 的 空间 分 布 
Fig.4 Decadal spatial distributions of longest duration over 3H Plain 
绝收 面积 数据 。 由 图 5 可 知 ,受灾 成 灾 和 绝收 面积 呈 波 动 下 降 趋势 。3 省 年 均 总 受灾 面积 由 1971 一 2000 年 
的 728.7 万 hm 下降 到 2000 一 2013 年 的 306.1 万 hm ,减少 了 58.0% ;成 灾 面 积 由 1971 一 2000 年 的 308.6 万 
hm2? 下 降 到 2000 一 2012 年 的 171.4 万 hm? ,减少 了 44.4% ;绝收 面积 由 1982 一 2000 年 的 64.2 万 hm 下 降 到 
2000 一 2012 年 的 32.6 万 hm ,减少 了 49.1%。 干旱 面积 显著 下 降 一 方面 与 农业 科技 进步 有 关 , 也 与 干旱 频率 
降低 .持续 性 减弱 有 关 。 
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表 3 黄 淮海 平原 1961 一 2014a 持续 时 间 最 长 的 干旱 事件 
Table 3 Top 10 drought events ranked by duration in 3H Plain from 1961 to 2014 


序号 站 点 发 生 时 段 持续 月 数 SPEI-3 均值 
Rank Station Persistent period Duration (month) Mean SPEI-3 

1 廊坊 1975.03 一 1976.01 11 —1.63 

2 塘 沾 1988.12 一 1989.10 11 -1.50 

3 天 津 1968.02 一 1968.10 9 -1.71 

4 沧 洲 1968.02 一 1968.10 9 —2.02 

5 黄骅 1965.06 一 1966.02 9 —1.65 

6 黄骅 1968.02 一 1968.10 9 -1.9 

7 南宫 1968.02 一 1968.10 9 —1.67 

8 石家庄 1968.02 一 1968.09 8 —1.44 

9 石家庄 1972.05 一 1972.12 8 —1.65 

10 Pb 台 1972.05 一 1972.12 8 -1.52 


350 


-= 加 


280 


210 


140 


70 


受灾 面积 
Drought area causing 10%/(X104hm2) 


成 灾 面积 
Drought area causing 30%/ (X 104hm2) 


1970 1976 1982 1988 1994 2000 2006 2012 1970 1976 1982 1988 1994 2000 2006 2012 


80 


60 
Hl 一 河南 
40 ! 9 ne 


20 


绝收 面积 
Drought area causing 80%/ (X 104hm2) 
3 ~---e>- 


1970 1976 1982 1988 1994 2000 2006 2012 


5 ” 河南、 河北 和 山东 省 的 农作物 干旱 受灾 面积 成 灾 面 积 和 绝收 面积 时 间 序 列 特征 


Fig.S Temporal evolution of drought area causing 10%, 30% and 80% yield loss in Henan, Hebei and Shandong province 


为 定量 分 析 农 业 实 际 干旱 面积 与 气象 干旱 指数 的 关系 ,计算 了 干旱 面积 与 SPEI 指数 的 Pearson 相关 系 
数 , 结 果 如 图 6。SPEI 指数 与 华北 三 省 受灾 .成 灾 绝收 面积 的 Pearson 相关 系数 7 为 负 值 ,表明 随 着 SPEI 指 
数 的 升 高 ( 趋 湿 ) ,人 遭受 旱灾 的 面积 (受灾 、 成 灾 和 绝收 面积 ) 减 少 ;从 7r 的 绝对 值 大 小 来 看 ,1r| 随 着 时 间 尺 度 和 
月 份 的 加 大 而 增 大 ,其 中 山东 省 12 月 份 12 个 月 时 间 尺 度 (Dec-SPEI- 12 ) 与 受灾 .成 灾 和 绝收 面积 的 171 分 别 
为 0.7.0.8 和 0.8 ,河南 省 Dec-SPEI- 12 与 受灾 .成 灾 和 绝收 面积 的 1r1 分 别 为 0.5 .0.6 和 0.6, 河 北 省 Dec-SPEI- 
12 与 受灾 .成 灾 和 绝收 面积 的 171 分 别 为 0.5.0.7 和 0.7, 这 表明 该 指数 能 够 反映 研究 区 域 的 干旱 与 旱灾 程度 。 
男 外 ,从 SPEI 指数 与 干旱 面积 相关 系数 可 以 看 出 ,相关 系数 + 在 7 月 份 以 后 趋 于 稳定 ,7 月 份 某 些 时 间 尺 度 的 
Irl 接 近 甚 至 大 于 Dec-SPEI- 12 所 对 应 的 1r| ,这 是 由 于 黄 淮海 地 区 降水 主要 集中 在 夏季 (7 一 9 月 )59 ,因此 
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7 月 份 之 前 的 干旱 灾害 损失 基本 上 决定 了 全 年 的 损失 。 仍 需 注意 的 是 ,受灾 、 成 灾 和 绝收 面积 分 别 代表 的 是 
因 灾 减产 10% 30% 和 80% 的 农作物 播种 面积 ,而 黄 淮 海地 区 旱灾 主要 发 生 在 冬小麦 生长 季 , 抽 穗 至 灌浆 
(5 一 6 月 ) 又 是 冬小麦 产量 形成 关键 期 ,因此 6 月 份 2 个 月 时 间 尺 度 SPEI( Jun-SPEI-2) 与 受灾 ,成 灾 和 绝 
收 面积 的 相关 系数 7 亦 比较 高 。 
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图 6 河南 、 河 北 和 山东 农作物 干旱 面积 与 不 同月 份 各 时 间 尺 度 的 SPEI 序列 的 Pearson 相关 系数 
Fig.6 Pearson correlation coefficients (+r) between monthly SPEI series at 1 一 12 month scales and drought area of Henan, Hebei and 


Shandong province 


3 讨论 

(1) 本文 改 用 基于 Penman-Monteith 蒸 散 模型 的 SPEI 指数 分 析 了 黄 淮 海平 原 近 52 年 来 气象 干旱 的 变化 
特征 ,指出 黄 淮 海平 原 1963 一 2014 年 干旱 指数 的 变化 趋势 呈现 明显 的 上 升 趋势 ,这 与 Kai Xu 等 人 !' 趾 基于 
Penman-Monteith 模型 的 SPEI 对 我 国 1961 一 2012 年 的 干旱 变化 趋势 的 研究 结果 一 致 ,而 与 其 他 学 者 "中 采 
用 原始 SPEI( 用 SPEI, 表示 , 蔡 散 模型 为 Thornthwaite) 指数 对 华北 地 区 的 干旱 趋势 结论 有 差异 。Meixiu Yu 等 
人 "中 对 全 国 609 个 气象 站 点 的 SPETI,, 分 析 表 明黄 淮海 地 区 特别 是 黄 淮 海北 部 呈现 显著 的 干旱 化 趋势 , 周 丹 
等 (由 用 同样 的 方法 分 析 得 出 华北 地 区 1962 一 2011 年 干旱 强度 加 重 、 干 旱 频 率 上 升 结论 , 李 伟 光 等 (利用 
SPEI%, 对 中 国 干旱 趋势 的 研究 亦 表明 我 国 普 遍 存 在 干旱 化 的 事实 ,西部 华北 和 东北 地 区 干旱 化 最 为 显著 。 
显然 ,不 同 的 蒸 散 模型 对 SPEI 指数 的 趋势 变化 产生 了 不 同 甚 至 反 转 的 效果 。 

研究 指出 ,Thornthwaite 与 Penman-Monteith 方法 计算 的 潜在 蒸 散 量 在 黄 淮海 地 区 的 历史 变化 呈现 出 了 两 
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种 相反 的 趋势 , 即 Thornthwaite 落 散 量 呈 显著 上 升 趋 势 ,而 Penman-Monteith 蒸 散 量 呈 显 著 下 降 趋势 “的 。 从 
SPEI 计算 过 程 中 可 以 看 出 ,干旱 指数 的 变化 趋势 由 降水 量 P 和 潜在 蒸 散 量 ET 共同 决 定 ,在 相对 干燥 的 黄 淮 
海地 区 ,潜在 蒸 散 量 RT, 通常 大 于 降水 量 P,SPEI 指数 的 变化 趋势 由 潜在 蒸 散 量 RT, 主导 ' 沾 。 因 此 ,尽管 该 地 
区 降水 量 呈 波动 下 降 '” ,基于 Penman-Monteith 蒸 散 模 型 的 SPEI 仍 呈 湿润 变化 。 男 外 ,有 报道 指出 采用 PDSI 
虽 标 和 标准 化 降水 指数 SPI 在 华北 地 区 的 干旱 趋势 是 呈 暧 干 化 趋势 的 '*”] ,这 主要 是 由 于 PDSI 干旱 指数 其 
水 分 平衡 中 莹 散 量 的 计算 采用 的 是 Thormmthwaite 模型 ,该 模型 计算 的 潜在 蒸 散 量 呈 上 升 趋势 导致 了 干旱 化 趋 
势 的 出 现 , 而 SPI 指数 仅 考虑 降水 的 影响 ,该 地 区 降水 多 年 气候 倾向 率 呈 下 降 趋 势 ,也 导致 了 干旱 化 的 趋势 。 
目前 ,对 基于 两 种 潜在 蒸 散 模型 的 SPEI 指数 在 我 国 的 对 比分 析 , 不 同学 者 得 出 较为 一 致 的 结论 , 即 
Penman-Monteith 模型 被 认为 是 最 具 物 理 意 义 和 最 可 靠 的 方法 ,并 且 经 常用 于 对 经 验 模型 有 效 性 的 验证 '*”3。 
研究 指出 ,由 于 Thornthwaite 模型 仅 考 虑 温度 变化 对 潜在 落 散 发 的 贡献 ,舍弃 了 对 藻 散 有 决定 性 作用 的 季节 性 
要 素 ,如 风速 .相对 温度 日照 时 数 等 ,该 模型 在 燕 散 量 较 低 的 地 区 高 舍 了 67, ,例如 在 我 国 东 南 地 区 ,而 在 蒸 
散 量 较 高 的 ,如 华北 地 区 则 低估 了 ET, ,无 法 有 效 捕 获 潜 在 蒸 散 的 空间 分 布 和 时 间 演 变 趋 势 , 其 在 中 国 地 区 的 
有 效 性 是 值得 怀疑 的 ' 站 。Penman-Monteith 公式 综合 考虑 了 热量 和 空气 动力 两 个 因子 ,在 估算 中 国 北方 的 潜 
在 蒸 散发 时 ,辐射 因子 项 和 空气 动力 因子 项 的 影响 都 不 能 忽视 ,并 且 随 着 气候 变 暖 ,其 中 空气 动力 因子 对 总 洪 
在 莹 散发 贡献 明显 增加 ,特别 是 在 北方 地 区 品 ] 。 刘 珂 和 姜 大 膀 ! 宫 对 比分 析 了 基于 两 种 潜在 蒸 散 发 算法 的 
SPEI 在 中 国 地 区 的 差异 ,指出 由 于 空气 动力 项 对 冬 .春季 北方 潜在 蒸 散发 的 贡献 显著 增加 , 基于 Penman- 
Monteith 公式 的 SPEI 相对 而 言 能 更 合理 地 描述 干 湿 变化 特征 。Kai Xu 等 人 :六 同样 指出 Thornthwaite 模型 过 


高 估计 了 气温 的 影响 。 

(2) 本 文 结合 历史 灾情 数据 ,首次 分 析 了 气象 干旱 指数 与 实际 农业 干旱 灾情 的 关系 。 黄 淮海 平原 降水 主 
要 集中 在 夏季 (7 一 9 月 ) ,空间 分 布 呈 南 多 北 少 的 格局 ,因此 该 地 区 以 春 旱 为 主 , 秋 早 频率 也 较 高 ,对 冬小麦 生 
= 影响 较为 显著 。 从 水 分 收 支 的 角度 来 看 ,农作物 的 实际 可 利用 水 量 为 有 效 降水 ( Pe) 与 实际 燕 散发 (ET) 差 
额 (Pe-57) 的 累计 值 , 因 此 基于 湾 在 燕 散 的 SPEI 指数 并 不 能 够 有 效 反映 实际 干旱 情况 。 因 此 , 想 要 进一步 识 
别 不 同 作物 生长 过 程 中 的 干旱 特征 ,应 从 农业 干旱 的 角度 入 手 ,在 计算 潜在 燕 散 量 时 订正 为 具体 作物 的 实际 
蔡 散 量 。 研 究 中 为 了 计算 作物 需 水 量 ,通常 采用 作物 系数 法 对 潜在 燕 散 量 进行 订正 "” 。 联 合 国 粮农 组 织 推 
荐 了 作物 系数 的 计算 方法 和 标准 状态 下 ( 白天 平均 最 低 相对 湿度 45% ,平均 风速 2 m/s, 半 湿润 气候 条 件 ) 各 
类 作物 的 作物 系数 参考 值 , 但 由 于 作物 系数 受 土壤 、 气 候 、 作 物 生 长 状况 和 栽培 管理 方式 等 诸多 因素 影响 ,对 
FAO 推荐 的 作物 系数 应 根据 当地 实际 情况 进行 订正 。 


4 结论 


(1) 蒸 散 模型 改 用 Penman-Monteith 公式 后 ,4 个 时 间 尺 度 (1.3.6 12 个 月 ) 的 SPEI 指数 在 黄 淮海 平原 均 
呈 上 升 趋势 , 即 与 采用 Thornthwaite 营 散 模型 的 SPEI 趋 于 干旱 化 不 同 ,新 的 干旱 指数 呈 湿 润 的 变化 趋势 。 这 
主要 是 由 于 两 种 蒸 散 模 型 计算 的 ET 趋势 差异 导致 的 ,基于 Penman-Monteith 公式 的 潜在 蒸 散 量 在 黄 淮海 平 
原 呈 显著 下 降 的 趋势 ,与 实际 蒸发 显 一 致 , 即 华 北 地 区 的 “蒸发 悖 论 " 现 象 。 

(2) 从 干旱 指数 的 历史 演变 来 看 ,基于 Penman-Monteith 蒸 散 公式 的 SPEI 指数 对 黄 淮海 平原 典型 干旱 年 
份 具 有 和 良好 的 适用 性 ; 黄 淮海 平原 6 个 农业 亚 区 的 干旱 频率 年 代 际 差异 表现 一 致 , 均 为 1960s 干旱 发 生 频 率 
最 高 ,而 2000 一 2014 年 最 低 ; 黄 淮海 平原 干旱 的 发 生 具 有 明显 的 持续 性 特点 ,1960s 干旱 持续 时 间 最 长 , 而 
2000 一 2014 年 为 历史 最 低 。 

(3) 本 文 利用 Pearson 相关 分 析 首 次 探讨 了 气象 干旱 指数 与 农业 实际 干旱 面积 的 关系 ,结果 表明 无 论 从 
相关 方向 还 是 相关 系数 的 大 小 ,SPEI 指数 对 农业 干旱 面积 的 变化 具有 良好 的 表征 能 力 , 月 值 干旱 指数 Dec- 
SPEI- 12 与 河南 .河北 和 山东 三 省 的 受灾 .成 灾 和 绝收 面积 相关 系数 绝对 值 最 大 。 
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